Conclusions

We have analyzed in detail the renormalization of the axial charge in nuclei by evaluating the matrix elements of the axial charge operator with relativistic wave functions, solutions of the Dirac equation with the relativistic potential, and with non relativistic wave function, solutions of the Schrödinger equation with an equivalent non relativistic potential. We have found renormalization effects due to the use of the relativistic wave functions, enhancing the axial charge in the direction found in earlier perturbative approaches for nuclear matter. However, the quantitative results differ from the estimates derived for nuclear matter significantly. Using the G-matrix derived from a realistic meson exchange model of the NN interaction [#rupr##1###] a perturbative estimate of the heavy meson exchange current contribution to the axial charge at nuclear matter saturation density [#7a##1###] would yield a renormalization factor of 1.4 and a non perturbative treatment would lead to enhancement as large as 1.7. For finite nuclei the enhancement factors considerably smaller, of the order of 1.2 - 1.3. We argue that this reduction of the renormalization effect is due to the smaller densities at the surface of finite nuclei, which are relevant for the evaluation of actual matrix elements. From these considerations we can also understand the dependence of the renormalization factor on the momentum transfer and on the transition actually considered.

The amount of axial charge renormalization depends on the model for the NN interaction. We have employed a relativistic meson exchange model (Potential version A of the Bonn potential [#rupr##1###]), which has been derived to reproduce NN scattering data. It is fair to quote at this point that using this potential in the present case there is the assumption that the relativistic potential constructed to reproduce NN scattering of on shell nucleons can be extrapolated to deal with negative energy states and on shell and off shell conditions. This is certainly a strong assumption from which all the microscopically constructed relativistic potentials suffer, and indeed different parametrizations of the NN amplitude on shell lead to different relativistic potentials [#11##1###]. Some efforts have been done to constrain the relativistic potential to be consistent with the #tex2html_wrap_inline1414#N elementary amplitudes [#12##1###] and this leads to potentials like the one obtained here but about one half their strength. Even then this potential is constructed at the level of the impulse approximation or low density limit, , and many body effects should modify it. It is clear that many efforts are still necessary to be able to claim that an unambiguous microscopical relativistic potential has been determined. On the other hand one can take a more phenomenological approach and say that a certain relativistic potential has a wide degree of phenomenological success, providing fair nuclear binding energies, spin-orbit splitting, nucleon nucleus cross sections and polarization observables, etc. [#10##1###,#klein##1###]. The potential we have used is one of such and provides empirical support for the axial charge renormalization found, but this does not exclude the possibility of other potentials with the same degree of phenomenological success and still providing different axial charge renormalization. The ultimate answer to this question is tied to the progress in our understanding of the meaning and accurate strength of the relativistic potential. Meanwhile, by using a fair and plausible model we have done detailed calculations and showed that the results are sufficiently different from the perturbative results to encourage the use of the present approach in future works dealing with the problem.



Two of us, A. Gil and E. Oset wish to acknowledge the hospitality of the University of Tübingen and H. Müther the one of the University of Valencia. E. Oset acknowledges support from the Humboldt Foundation. The work has been partially supported by the EU, program, Human Capital and Mobility contract no. CHRX-CT 93-0323, the CICYT contract no. AEN 93-1205 and the Graduiertenkolleg ``Struktur und Wechselwirkung von Hadronen und Kernen'' of the Deutsche Forschungsgemeinschaft (DFG Mu 705/3)

<#370#>Appendix: Matrix elements of the axial charge operator.<#370#>


<#371#>A) Relativistic case:<#371#> We write here the matrix element for the #math57#γ0γs operator between relativistic wave functions

#math58#

;SPMlt; n'l'j'm'| γ0γ5ei#tex2html_wrap_indisplay1418##tex2html_wrap_indisplay1419#| nljm ;SPMgt; #tex2html_wrap_indisplay1420#(a.1)

We distinguish 4 cases

#math59#

a)j' = l' + 1/2, j = l + 1/2
b)j' = l' + 1/2, j = l - 1/2
c)j' = l - 1/2, j = l + 1/2
d )j' = l - 1/2, j = l - 1/2

and the resulting matrix element is

#math60#

#tex2html_wrap_indisplay1431##tex2html_wrap_indisplay1432#iλ(- i)#tex2html_wrap_indisplay1433#r2dr#tex2html_wrap_indisplay1434##tex2html_wrap_indisplay1435##tex2html_wrap_indisplay1436# - #tex2html_wrap_indisplay1437##tex2html_wrap_indisplay1438##tex2html_wrap_indisplay1439#jλ(qr)

#math61#

(2λ +1)1/2Yλ, m'-m*(#tex2html_wrap_indisplay1441#)Ai#tex2html_wrap_indisplay1442#(a.2)

where Ai is given for each of the cases a) b) c) d) listed above by

#math62#

Aa = #tex2html_wrap_indisplay1445#C(j + 1/2, λ, j' - 1/2;000)

#math63#

{(j' + m')1/2(j + 1 - m)1/2C(j + 1/2, λ, j' - 1/2;m - 1/2, m' - m)

#math64#

- (j' - m')1/2(j + 1 + m)1/2C(j + 1/2, λ, j' - 1/2;m + 1/2, m' - m)}#tex2html_wrap_indisplay1448#(a.3)

#math65#

Ab = #tex2html_wrap_indisplay1450#C(j - 1/2, λ, j' - 1/2;000)

#math66#

{(j' + m')1/2(j + m)1/2C(j - 1/2, λ, j' - 1/2;m - 1/2, m' - m)

#math67#

+ (j' - m')1/2(j - m)1/2C(j - 1/2, λ, j' - 1/2;m + 1/2, m' - m)#tex2html_wrap_indisplay1453#(a.4)

#math68#

Ac = #tex2html_wrap_indisplay1455#C(j + 1/2, λ, j' + 1/2;000)

#math69#

{(j' + 1 - m')1/2(j + 1 - m)1/2C(j + 1/2, λ, j' + 1/2;m - 1/2, m' - m)

#math70#

+ (j' + 1 + m')1/2(j + 1 + m)1/2C(j + 1/2, λ, j' + 1/2;m + 1/2, m' - m)}#tex2html_wrap_indisplay1458#(a.5)

#math71#

Ad = #tex2html_wrap_indisplay1460#C(j - 1/2, λ, j' + 1/2;000)

#math72#

{(j' + 1 - m')1/2(j + m)1/2C(j - 1/2, λ, j' + 1/2;m - 1/2, m' - m)

#math73#

- (j' + 1 + m')1/2(j - m)1/2C(j - 1/2, λ, j' + 1/2;m + 1/2, m' - m)}#tex2html_wrap_indisplay1463#(a.6)


<#434#>B) Non relativistic case:<#434#> we evaluate matrix elements of the #math74##tex2html_wrap_inline1465#(#tex2html_wrap_inline1466# + #tex2html_wrap_inline1467#)/2M operator between non relativistic states

#math75#

;SPMlt; n'l'j'm'|#tex2html_wrap_indisplay1469#(#tex2html_wrap_indisplay1470# + #tex2html_wrap_indisplay1471# ')/2M| nljm ;SPMgt; #tex2html_wrap_indisplay1472#(a.7)

The derivation of this matrix elements requires a bit more algebra than the non relativistic case. With the help of some useful formulas from the appendix of ref. [#13##1###] we obtain the following result

#math76#

;SPMlt; n'j'l'm'|#tex2html_wrap_indisplay1474##tex2html_wrap_indisplay1475#ei#tex2html_wrap_indisplay1476# #tex2html_wrap_indisplay1477#| nljm ;SPMgt; =

#math77#

i(- 1)j'+l'+1/2#tex2html_wrap_indisplay1479##tex2html_wrap_indisplay1480##tex2html_wrap_indisplay1481#C(jλj';1/2, 0, 1/2)

#math78#

C(jλj';m, m' - m)iλYλ, m'-m*(#tex2html_wrap_indisplay1483#)#tex2html_wrap_indisplay1484##tex2html_wrap_indisplay1485##tex2html_wrap_indisplay1486#

#math79#

F(n'l'j', nlj;λ)#tex2html_wrap_indisplay1488#(a.8)

with #math80#λ + l + l' an odd number, where the last function is given by

#math81#

F(n'l'j', nlj;λ) =

#math82#

δj, l+1/2(l + 1)1/2#tex2html_wrap_indisplay1492#r2drφl'(r)[#tex2html_wrap_indisplay1493# - #tex2html_wrap_indisplay1494#φl(r)]jλ(qr)

#math83#

- (- 1)j-j'+l-l'δj', l'+1/2(l' + 1)1/2#tex2html_wrap_indisplay1496##tex2html_wrap_indisplay1497##tex2html_wrap_indisplay1498#

#math84#

#tex2html_wrap_indisplay1500#r2dr[#tex2html_wrap_indisplay1501# - #tex2html_wrap_indisplay1502#φl'(r)]φl(r)jλ(qr)

#math85#

- δj, l-1/2l1/2#tex2html_wrap_indisplay1504#r2drφl'(r)[#tex2html_wrap_indisplay1505# + #tex2html_wrap_indisplay1506#φl(r)]jλ(qr)

#math86#

+ (- 1)j-j'+l-l'δj', l'-1/2l'1/2#tex2html_wrap_indisplay1508##tex2html_wrap_indisplay1509##tex2html_wrap_indisplay1510#

#math87#

#tex2html_wrap_indisplay1512#r2dr[#tex2html_wrap_indisplay1513# + #tex2html_wrap_indisplay1514#φl'(r)]φl(r)jλ(qr)#tex2html_wrap_indisplay1515#(a.9)